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Brownian dynamics simulations of a coarse-grained bead-spring chain model, with Debye-Hiickel electrostatic 
interactions between the beads, are used to determine the root-mean-square end-to-end vector, the radius 
of gyration, and various shape functions (defined in terms of eigenvalues of the radius of gyration tensor) 
of a weakly-charged polyelectrolyte chain in solution, in the limit of low polymer concentration. The long¬ 
time diffusivity is calculated from the mean square displacement of the centre of mass of the chain, with 
hydrodynamic interactions taken into account through the incorporation of the Rotne-Prager-Yamakawa 
tensor. Simulation results are interpreted in the light of the OSFKK blob scaling theory (R. Everaers, 
A. Milchev, and V. Yamakov, Eur. Phys. J. E 8, 3 (2002)) which predicts that all solution properties 
are determined by just two scaling variables—the number of electrostatic blobs X, and the reduced Debye 
screening length, Y. We identify three broad regimes, the ideal chain regime at small values of Y, the blob- 
pole regime at large values of Y, and the crossover regime at intermediate values of Y, within which the 
mean size, shape, and diffusivity exhibit characteristic behaviours. In particular, when simulation results are 
recast in terms of blob scaling variables, universal behaviour independent of the choice of bead-spring chain 
parameters, and the number of blobs X, is observed in the ideal chain regime and in much of the crossover 
regime, while the existence of logarithmic corrections to scaling in the blob-pole regime leads to non-universal 
behaviour. 
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I. INTRODUCTION 

Many synthetic polymers and most biopolymers are 
polyelectrolytes. Their use in a range of industrial and 
biological applications makes a thorough understanding 
of their behaviour highly desirable. While the scaling of 
many important observable properties of neutral poly¬ 
mers in solution has been successfully predicted by blob 
theorieSfii^ much is yet to be understood with this ap¬ 
proach even for the simple case of dilute polyelectrolyte 
solutions at equilibrium. The electrostatic blob, which 
sets the length scale at which the energy of electrostatic 
interactions is of order ksT (where ks is the Boltzmann 
constant and T is the temperature), provides the ba¬ 
sis for scaling theories of such a system^ A commonly 
used scaling theory for polyelectrolyte solutions is the 
OSFKK scaling picture^^— (after Odjik, Skolnick, Fix- 
man, Khokhlov, and Khachaturian), which predicts the 
scaling of the mean chain size in different scaling regimes 
based on the dominating physics. Even though simula¬ 
tions have examined the predictions of OSFKK scaling 
theory in detail;^ many questions remain unanswered and 
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not many studies have applied the blob picture to proper¬ 
ties other than the root-mean-square end-to-end vector. 
The aim of this work is to carry out Brownian dynam¬ 
ics (BD) simulations in order to investigate the depen¬ 
dence of several static and dynamic properties of dilute 
polyelectrolyte solutions on the intrinsic parameters that 
govern their behaviour. In particular, in addition to the 
mean size, we use the OSFKK scaling picture to inter¬ 
pret simulation predictions of the shape and diffusivity of 
polyelectrolyte chains. A variety of shape functions, de¬ 
fined in terms of the eigenvalues of the radius of gyration 
tensor, are used to examine changes in polymer shape in 
the different regimes of conformational phase space. Fur¬ 
ther, hydrodynamic interactions have been incorporated 
via the Rotne-Prager-Yamakawa tensor ,i2i£ in order to ob¬ 
tain an accurate prediction of the chain diffusivity in var¬ 
ious scaling regimes. The OSFKK scaling picture ignores 
the presence of logarithmic corrections in the blob-pole 
regime of the conformational phase diagram, which are 
predicted to be important by more refined scaling theo¬ 
ries'^ For all the properties examined here, we explore 
whether logarithmic corrections can account for depar¬ 
tures from OSFKK scaling in the blob-pole regime. 

In most scaling theories, polyelectrolyte molecules are 
represented as freely-jointed chains with Npi Kuhn steps, 
each of length 6 k- The number of charges per Kuhn step. 
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FIG. 1. The phase diagram in the OSFKK scaling pictnre.— 
Various regimes are marked I-VI, with scaling predictions for 
the end-to-end vector re, in each regime, listed in Table [H 
Equations governing the boundaries between regimes are in¬ 
dicated in the figure, and the abbreviation “RW” implies ran¬ 
dom walk statistics are obeyed. 


/ = aA^m.K, where, A^m.K is the number of monomers per 
Kuhn step, and a is the degree of ionization per chain, 
is used to characterize the extent of ionic group dissocia¬ 
tion. Counterions and salt ions are not modelled explic¬ 
itly, and the screening of monomer charges due to the 
presence of free ions is quantified via the Debye screen¬ 
ing length, /d- Another important length scale is the 
Bjerrum length. Is, which is the distance at which the 
Coulomb energy between two unit elementary charges in 
the solvent is equal to the thermal energy ksT, and is 
defined by, 

^ 47reo£rfcB2^ ^ ^ 

where e is the charge of an electron, eg is the vacuum per¬ 
mittivity, and £r is the relative dielectric constant of the 
solvent. Such a representation of a dilute polyelectrolyte 
solution is referred to here as the bare-model, and in what 
follows, it is always assumed that the macromolecules are 
dissolved in a theta solvent (which is simulated by switch¬ 
ing off excluded volume interactions). 

The electrostatic blob denotes the length scale below 
which the conformation of a polyelectrolyte chain is prac¬ 
tically unaffected by electrostatic interactions, since the 
electrostatic energy of the sub-chain within a blob is less 
than the thermal energy. This length scale can be used 
to divide the chain into a number blobs, X, of diam¬ 
eter, ^ei, with random walk statistics followed within 
the blob, while the conformation of the chain of blobs 
as a whole depends on the electrostatic interactions be¬ 
tween the blobs. In addition, within the OSFKK scaling 
picture, the screening of electrostatic interactions is ac¬ 
counted for by the scaled variable V, which is the ratio 
of the Debye length to the blob size. In terms of the 
bare-model parameters, the OSFKK expressions for the 


blob scaling variables are^. 


and, 

Y=^ (3) 

Cel 

where, Id = and the scaled blob size ^ei = Cei/^K, 

is given by, 

( 4 ) 

with, Zb = In/bK- Note, for future reference, that it is 
sufficient to prescribe the reduced set of bare-model pa¬ 
rameters {Nk, JjIbjId} in order to determine the scal¬ 
ing variables X and Y. The OSFKK scaling picture is 
briefly summarized in Fig. [U which is a phase diagram in 
(A, y) space, highlighting the different scaling regimes. 
Table U lists the scaling expressions for the chain size, in 
terms of both blob and bare parameters, in each of the 
regimes. In regime I, when A <C I (the blob is much 
larger than the chain), re/^ei oc A^^^, where re is the 
root-mean-square end-to-end vector. The Debye length 
is irrelevant in this regime as thermal energy completely 
dominates the electrostatic interactions between the seg¬ 
ments. When A ^ 1 (when the chain is composed of 
many electrostatic blobs) the XY plane is divided into 
different regimes according to how large Y is compared 
to A. The blob-pole conformation occurs in regime II, 
when Y —> oo, and electrostatic interactions between 
blobs are unscreened. In this regime, re/^ei oc A. When 

Y decreases below A^^^, electrostatic interactions act as 
short ranged repulsive forces and a crossover regime fol¬ 
lows, which, according to Khokhlov and Khachaturian^ 
is decomposed into regimes III-V. Here the confirmations 
are governed by an intricate interplay of chain stiffness 
and excluded volume interactions (both of electrostatic 
origin). Since our numerical investigations lack the res¬ 
olution to study the details of these regimes we collec¬ 
tively refer to them as the crossover regime. Finally when 

Y <C the Debye length is sufficiently small for the 

scaling to return to that of an ideal chain. 

Everaers, Milchev, and Yamakov— have carried out ex¬ 
tensive single chain Monte Carlo simulations with Debye- 
Hiickel electrostatic interactions and their results support 
the OSFKK scaling picture for the end-to-end vector. 
Pattanayek and Prakash ti used single chain Brownian 
dynamics simulations with Debye-Hiickel electrostatic in¬ 
teractions in order to investigate the scaling of the end- 
to-end vector and viscometric functions in simple shear 
flow. They showed that when Brownian dynamics sim¬ 
ulation data were represented in terms of blob scaling 
variables, predictions independent of the choice of pa¬ 
rameters in the bead-spring chain model were obtained. 
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TABLE 1. Scaling expressions for the end-to-end vector in the various regimes shown schematically in Fig. [T1 in terms of bare 
and blob parameters (the value of the Flory excluded volume exponent has been approximated to be 3/5). The variables X 
and Y are defined in Eqs. © and respectively. 


Regime 

re/Cel (blob parameters) 

relbK (bare parameters) 

I 

Xi/2 


II 

X 

j2/3 fl/3 

III 

YX^t2 

/A3 fi/3 Ari/2 

IV 

y^3/5j^3/5 

.8/15 74/15 p/5 .p/5 

V 

y2/5j^3/5 


VI 

yl /2 



As in the present study, these Brownian dynamics simu¬ 
lations were unable to distinguish between the different 
regimes in the crossover region due to the computational 
cost of simulating the extremely long chains required to 
verify OSFKK scaling. By including a characteristic non- 
dimensional shear rate as an additional scaling variable 
in shear flow, Pattanayek and Prakash— were able to 
show that the equilibrium electrostatic blob model pro¬ 
vides a framework to obtain universal scaling even for 
non-equilibrium properties. Both these simulation stud¬ 
ies, however, did not attempt to describe the scaling of 
other static properties, or the diffusivity in the different 
regimes, nor did they examine the occurrence of logarith¬ 
mic corrections in the blob-pole regime. 

The appropriateness of using a Debye-Hiickel poten¬ 
tial to represent electrostatic interactions in Brown¬ 
ian dynamics simulations has been examined by Stoltz, 
de Pablo, and Grahami^ They compared the predic¬ 
tions of chain size as a function of Bjerrum length (at 
various concentrations in the dilute regime) of a bead¬ 
spring chain model with explicit counterions and pair¬ 
wise Coulomb interactions between charges, with those 
of a model with pair-wise Debye-Hiickel electrostatic in¬ 
teractions between beads on chains. They found that the 
results of both models are nearly identical for values of 
Bjerrum length roughly equal to or less than the distance 
between the beads on a chain (all of which were assumed 
charged). For larger values of Bjerrum length, in ac¬ 
cord with Manning’s theorythey observed the onset 
of counterion condensation. We have adopted a Debye- 
Hiickel potential in this work, and as discussed in greater 
detail subsequently, chosen parameter values in the bead¬ 
spring chain model to ensure that the simulations always 
remain in the regime where this approximation is valid, 
and there is no counterion condensation. 

Property predictions from Brownian dynamics simu¬ 
lations of bead-spring chains depend on several model 
parameters, such as the number of beads, the charge 
on the beads, the finite extensibility parameter for the 
springs, and so on. A key issue is a rational choice of val¬ 
ues for these parameters. The earlier BD simulations of 
mean size and viscometric functions by Pattanayek and 
Prakash— has shown the advantage of using blob scal¬ 


ing variables to interpret results of simulations since this 
leads to a description independent of the level of coarse- 
graining, which is extremely useful for comparing simu¬ 
lation results with experiments. In this study, we inves¬ 
tigate the shape and diffusivity of polyelectrolyte chains 
in the various regimes of the phase diagram, by recasting 
BD simulation results in terms of blob scaling variables, 
and examine their independence from the specific choice 
of bead-spring chain model parameters. 

Many studies have shown that the shape of a neutral 
polymer chain, even at equilibrium, is not spherical about 
the centre of mass of the chainji^^— The nature of the 
asymmetry in chain shape has been examined in terms 
of a number of different quantities, such as the degree of 
prolateness, the asphericity, the acylindricity, the shape 
anisotropy, and so on, which are functions of the eigen¬ 
values of the radius of gyration tensor, since the breaking 
of symmetry is reflected in the three eigenvalues differing 
from each other. The symmetry or otherwise of polyelec¬ 
trolyte chain shapes, particularly in the different scal¬ 
ing regimes, has not yet been systematically investigated. 
Here, we examine if a universal description of polyelec¬ 
trolyte chain shapes can be obtained, when BD simula¬ 
tion results are represented in terms of electrostatic blob 
scaling variables. 

While the OSFKK scaling picture ignores the presence 
of logarithmic corrections to the scaling of chain size 
with degree of polymerisation in the blob-pole regime, 
their existence has been derived in a number of differ¬ 
ent ways, ranging from Flory type energy minimisation 
arguments,— to refined scaling theories that account for 
the nonuniform stretching of polyelectrolyte chains along 
the elongation axisiii^ In Appendix (for the sake of 
completeness), we have used a Flory type argument to 
show how logarithmic corrections arise in regime H. In 
particular, it can be shown that re obeys the following 
scaling expression in terms of the blob scaling variable 
A, 

g^A[lnA]^ (5) 

Liao, Dobrynin, and Rubinstein— have carried out 
molecular dynamic simulations of bead-spring chains 
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with explicit counterions, and have shown in terms of 
bare-model parameters, that for sufficiently long chains, 
the scaling of the end-to-end vector does indeed exhibit 
logarithmic corrections in the blob-pole regime. Here, we 
examine whether results of BD simulations in the blob- 
pole regime exhibit logarithmic corrections as described 
by Eq. ([S|). 

A common assumption in the various theoretical de¬ 
scriptions of the blob-pole regime is that electrostatic in¬ 
teractions lead to chain stretching along one direction, 
while leaving the chain conformation unperturbed in di¬ 
rections perpendicular to the stretching direction»iS In 
terms of blob scaling variables, this implies that chain di¬ 
mensions lateral to the stretching direction are expected 
to scale as ^ei A' 2 . By examining the eigenvalues of the 
gyration tensor, we verify if chain dimensions perpendic¬ 
ular to the stretching direction indeed obey ideal chain 
scaling laws. 

The concept of the Zimm diffusivity of a blob has been 
successfully used to develop scaling relations for the diffu¬ 
sivity of neutral polymer chains, both in the dilute con¬ 
centration regime (in terms of thermal blobs^), and in 
the semidilute regime (in terms of correlation blobs^i^). 
Here, we examine if the scaling of the diffusivity of poly¬ 
electrolyte chains in the various regimes of the phase dia¬ 
gram, becomes independent of bead-spring chain param¬ 
eters, when the Zimm diffusivity of an electrostatic blob 
is used to interpret simulation results. In particular, since 
the blob-pole is reminiscent of the shish-kebab model for 
rodlike polymers^^ (with blobs taking the place of beads), 
we examine if the diffusivity of a polyelectrolyte chain in 
regime B can be understood in terms of the translational 
diffusivity of rodlike polymers. 

The paper is structured as follows. In section B we 
describe the bead-spring chain model used to represent 
polyelectrolyte chains, and the governing equations for 
the time evolution of the position vectors for the beads. 
Section IB, which summarises our results and the rele¬ 
vant discussions, is subdivided into four sections; the first 
looks at the scaling of the chain size and the extent of log¬ 
arithmic corrections in the blob-pole regime, the second 
and third consider the scaling of various functions that 
describe the shape of the chain, while the fourth exam¬ 
ines the scaling of chain diffusivity and relaxation time. 
Finally, the key findings of this work are summarised in 
section IV. 


II. THE BEAD-SPRING CHAIN MODEL 

A. Governing equations and Brownian dynamics 
simulations 

A polyelectrolyte chain is modelled in the BD sim¬ 
ulations by a coarse-grained version of the bare-model, 
i.e., by a bead-spring chain consisting of beads of ra¬ 
dius a, connected linearly by (Vb — 1) finitely extensible 
non-linear elastic (FENE) springs, with spring constant 


H and maximum stretch Qq. The beads act as centres 
of frictional resistance, with a Stokes friction coefficient, 
C = Girris a (where rjs is the solvent viscosity). The total 
charge on the bead-spring chain is set equal to that of a 
chain in the bare-model, distributed uniformly along the 
length of the chain with each bead having an identical 
charge q, given by 


The time evolution of the position vector r^{t) of bead 
fi, is described by the non-dimensional stochastic differ¬ 
ential equation^ 

= ^ (7) 

where, the length scale Ih = \/k^TjH and time scale 
\h = C,I^H have been used for non-dimensionalization. 
The dimensionless diffusion tensor is a 3 x 3 ma¬ 
trix for a fixed pair of beads /i and v. It is related to 
the hydrodynamic interaction tensor, as discussed fur¬ 
ther subsequently. The sum of all the non-hydrodynamic 
forces on bead v due to all the other beads is represented 
by F,y, the quantity Wj, is a Wiener process, and 
is a non-dimensional tensor whose presence leads to mul¬ 
tiplicative noise Its evaluation requires the decompo¬ 
sition of the diffusion tensor. Defining the matrices X> 
and 3 as block matrices consisting ot N x N blocks each 
having dimensions of 3 x 3, with the (/i, z 2 )-th block of T> 
containing the components of the diffusion tensor 
and the corresponding block of B being equal to B^j,, 
the decomposition rule for obtaining B can be expressed 
as 

B B^ ( 8 ) 

The non-hydrodynamic forces on a bead /r are com¬ 
prised of the non-dimensional spring forces and 
non-dimensional electrostatic interaction forces F®®, i.e., 
F^ = F^ -I- F®®. The entropic spring force on bead 
/i due to adjacent beads can be expressed as F^ = 
F®(Q^) —F®(Q^_i) where F®(Q^_j) is the force between 
the beads ^ — 1 and /i, acting in the direction of the con¬ 
nector vector between the two beads Q^_i = — r^_i. 

Specifically, the spring force in the FENE springs used 
here is given by, 

where b = HQ^/k-QT is the dimensionless finite exten¬ 
sibility parameter. The vector F®® is given in terms of 
the dimensionless electrostatic potential C/®® between the 
beads /i and v of the chain. 
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The OSFKK scaling theory assumes Debye-Hiickel elec¬ 
trostatic interactions^^ between charges, which is justi¬ 
fied for weakly charged polyelectrolyte chains in the ab¬ 
sence of Manning counterion condensation! 37 

adopt a Debye-Hiickel potential in this work, 

= —“p(-T?) ( 11 ) 

with = |r^i,|, where is the vector 

between beads v and /r, and Zg and Zg are the nondimen- 
sional Bjerrum and Debye lengths, respectively. 

The non-dimensional diffusion tensor is related to 
the non-dimensional hydrodynamic interaction tensor fi 
through 

+ {1 - f2{r^ - r^) ( 12 ) 

where S and represent a unit tensor and a Kronecker 
delta, respectively, while f2 represents the effect of the 
motion of a bead /i on another bead v through the dis¬ 
turbances carried by the surrounding fluid. The hydro- 
dynamic interaction tensor fi is assumed to be given by 
the Rotne-Prager-Yamakawa (RPY) regularisation of the 
Oseen function 

rr 

17(r) = S f22~2 (1^) 

where for r = |r| > 2a*, 


3 a* 

( 2 a*^\ 

! ^ 3 a*, 

f a*2\ 


1 + 77 ^ 

and 122 = - — 

1-2^ 

4 r 

V 3 r^ J 

4 r ' 

i ^2 1 


(14) 

while for 0 < r < 2a*, 

= and 122 = ^ ^ (15) 

32 a* 32 a* 

Here, a* is the non-dimensional bead radius, which is 
related to the conventionally defined^ hydrodynamic in¬ 
teraction parameter, h*, by a* = y/nh*. We have set 
a* — 0.5 in all simulations in which hydrodynamic in¬ 
teractions are incorporated. A choice of h* close to 0.25 
ensures that the non-draining limit is reached (and con¬ 
sequently universal predictions obtained), at relatively 
small values of 

The spatial configuration of the chain at any time i, 
i.e., r^(t) for all beads fi = 1,..., N^, is obtained by inte¬ 
grating Eq. 0 using a semi-implicit predictor-corrector 
scheme proposed by Prabhakar and Prakash.lS In the 
presence of fluctuating HI, the problem of the computa¬ 
tional intensity of calculating the Brownian term is re¬ 
duced by the use of a Chebyshev polynomial represen¬ 
tation for the Brownian term4^i^ We have adopted this 
strategy, and the details of the exact algorithm followed 
here are given in Ref. 

While the majority of results reported here have been 
obtained with the “single-chain” BD algorithm described 
in Ref. IH, the predictions of eigenvalues and shape 
functions have been obtained with the “multi-chain” 


BD algorithm described in Ref. The most sig¬ 

nificant difference between the two algorithms is that 
the latter (in which multiple chains are simulated in a 
box with periodic boundary conditions), accounts for 
inter-particle hydrodynamic and electrostatic interac¬ 
tions in addition to intra-particle interactions. Essen¬ 
tially, Nc bead-spring chains are simulated in a cubic 
simulation box of length L, such that the concentration 
c = N\yNc/L^ ■ With the overlap concentration c* defined 
by, c* = Afb/[(4/3)7rr3 gq], where rg,eq is the radius of gy¬ 
ration of a chain at equilibrium, we maintain a scaled 
concentration c/c* = 10“®, to ensure that the system is 
in the dilute limit. Since the typical distance between 
chains at these concentrations is much greater than the 
Debye screening lengths considered here, the strength of 
short-ranged inter-molecular Debye-Hiickel electrostatic 
interactions is effectively zero. The equivalence of the 
use of either of the algorithms has been verified by com¬ 
parison of predictions of the end-to-end vector and the 
radius of gyration, and ensuring that no difference was 
observed. Plots which contain data from the multi-chain 
BD algorithm are identified as such in the figure captions, 
and unless explicitly stated, most plots report data ob¬ 
tained using the single chain BD algorithm. 

B. Size, shape and diflusivity 

The two static properties examined here are: (i) the 
end-to-end distance, r^ = \/ (rg), with, 

(^5 = ((rNb - ri) • (rNb - d)) (16) 

where, ((.)) represents an ensemble average, and and 
ri are the positions of the two beads at either end of 
the chain, and, (ii) the radius of gyration of the chain, 

»'g = with, 

(^g)=(A?) + (A^) + (A2) (17) 

where, A^, A^, and A§ are eigenvalues of the gyration 
tensor G (arranged in ascending order), with, 

^ JVb JVb 

^ ~ 22 ^ ^ ^ ( 1 ^) 
b ^=1 

Note that, G, Af, A|, and A§ are calculated for each 
trajectory in the simulation before the ensemble averages 
are evaluated. 

The asymmetry in equilibrium chain shape has been 
studied previously in terms of various functions defined 
in terms of the eigenvalues of the gyration tensor.— 
Apart from Xf, A 2 , and A§, themselves, we have exam¬ 
ined the following shape functions: the asphericity (B), 
the acylindricity (C), the degree of prolateness {S and 
5*), and the shape anisotropy (k^ and ), as defined 
in Table im As is evident from the Table, the latter two 
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TABLE II. Definitions of shape functions in terms of eigenvalues of the gyration tensor, G. Note that, 7i = Ai + A| + A 3 , and, 
I 2 — A 1 A 2 “t“ A 2 A 3 T A 3 A]^, are invariants of G. 


Shape function 

Definition 



Asphericity22 

B = (A 3 ) — - [(Af) - 1 - (A 2 )] 

(19) 


AcylindridtySS 

C = (A^) - (A?) 

( 20 ) 



/(3A?-7?)(3Ai-/?)(3Ai 

-ii)\ 

( 21 ) 


\ (/.")• 

/ 

Degree of prolateness^*^ 

^ ((3Af-7?)(3Ai-7?)(3Ai- 

-I?)) 

( 22 ) 


((/?)^) 



.--,-3(1) 

(23) 




(24) 



quantities are evaluated in terms of two different defi¬ 
nitions. Typically, it is easier to evaluate S and for 
analytical calculations, rather than S* and , which 
require the evaluation of averages of ratios of fluctuating 
quantities.— 

The only dynamic property that is directly calculated 
here is the long-time self-diffusion coefficient D, defined 
by, 


D = lim 

t—¥00 


|re(t)-re(0)p \ 
6t / 


(25) 


where. 


is the position vector of the 


= 

center of mass of the chain. 

The calculation of the radius of gyration and the dif- 
fusivity of a chain, enables an estimation of the longest 
relaxation time from the expression. 


(rl) 

D 


(26) 


An ensemble averaging method has been used to es¬ 
timate all the properties, details of which are given in 
Appendix A of Ref. H^. We find that trajectories level 
off after roughly one relaxation time, but a true station¬ 
ary state from which equilibrium properties can be esti¬ 
mated, generally occurs after tens of relaxation times. 
Typically, more than 0(10®) data points are used in 
the estimation of averages and standard errors of mean, 
where a data point is a saved property value at some time 


after the beginning of the stationary state. For proper¬ 
ties evaluated with the single chain code, a time step size 
of 0.005 was used, and data was saved roughly after every 
200 time steps. Averages were carried out over 1000 tra¬ 
jectories, with 1000 data points in each trajectory. For 
properties evaluated with the multi-chain code, a time 
step size of 0.005 was used, and data was saved roughly 
after every 10 to 20 time steps. Averages were carried out 
over 64 trajectories, each with 15 chains in the simulation 
box, and 2000 data points in each trajectory. 


C. Mapping between bead-spring chain and blob variables 

The concept of an electrostatic blob enables the rep¬ 
resentation of experimental and simulation data for 
dilute polyelectrolyte solutions in terms of a signifi¬ 
cantly reduced number of scaling variables ({X, F}), 
when compared to the number of bare-model parame¬ 
ters {{Nk, J/Ib/Id})- Further, the universal nature of 
property predictions when expressed in terms of blob 
scaling variables has been established previously by sim¬ 
ulations, which have shown that various combinations of 
bare-model parameters lead to identical results, provided 
values of the scaling variables are kept fixed In order 
to achieve a similar demonstration for BD simulations 
of bead-spring chains, it is necessary to map the bead¬ 
spring chain parameters {W, 5, q,/g,/p} onto the blob 
scaling variables {X,Y}. This issue has been discussed 
previously by Pattanayek and Prakashfli For the sake of 
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TABLE III. Values for the bead-spring chain parameters {Ab,6, the resultant scaling variables X and and the 

corresponding values for the bare model parameters {Vk, /, ^b}, used in Brownian dynamics simulations. Note that the listed 
values of the scaling variable Y are specified independently of other parameters by suitably choosing (not tabulated here). 
L* is the non-dimensional contour length of the bead-spring chain, and 70 is the Manning parameter. 


Blob model 

Bead-spring chain model 

Bare-model 

A 

Y 

a 

{Nb,6,g,ZS} 

L* 

{NkJM 

70 

20.3 

30 

2.07 

{32,72.42,1.1,0.69} 

272 

(800, 0.044, 2.09} 

0.092 

24 

0.01,0.05,0.1,0.5,1,1.6, 

3,4.8,10,20,30,100 

1.77 

(28,61.67,1.79,0.41} 

220 

(600,0.0833,1.15} 

0.096 

1.90 

(32,72.42,1.25,0.69} 

272 

(800,0.05,2.08} 

0.104 

2.21 

(40,61.67,1.38,0.4} 

314 

(867,0.0638,1.13} 

0.072 

27.6 

30 

1.77 

{32,61.67,1.79,0.41} 

251 

{689,0.0831,1.16} 

0.096 

30.74 

1.68 

(32,72.42,1.5,0.69} 

272 

(800,0.06,2.09} 

0.182 

34.83 

1.72 

(38, 72.42,1.45,0.69} 

323 

(955,0.0577,2.09} 

0.170 

37.67 

1.70 

(40, 72.42,1.48,0.69} 

340 

(1006,0.0588,2.09} 

0.125 

48 

0.01,0.05,0.1,0.5,1,1.6, 

3,4.8,10,20,30,100 

1.73 

{52,72,1.81,0.45} 

441 

{1310,0.0718,1.36} 

0.098 

1.85 

(60,61.67,1.70,0.41} 

471 

(1310,0.0777,1.16} 

0.090 

72 

0.01,0.05,0.1,0.5,1,1.6, 

3,4.8,10,20,30,100 

1.70 

{75,72,1.95,0.41} 

636 

{1900,0.0771,1.24} 

0.096 

1.80 

(85,61.67,1.61,0.5} 

668 

(1900,0.0732,1.41} 

0.103 


completeness, we summarise their arguments here. 

The key first step is to map bead-spring chain pa¬ 
rameters onto bare-model parameters, after which the 
mapping onto blob scaling variables is straight-forward. 
This is achieved by equating the end-to-end vectors and 
the contour lengths of the bare-model and bead-spring 
chains. If L is the contour length, then, 


L = Nk 6k = (-^b — 1) <3o = (M) — 1) (27) 

Further, using the analytical result for the end-to-end 
vector of FENE chains,— we get, 

(r2)=lVK6^ = (A^b-l)^;& (28) 

Equations (l27)) and (l28l) can be solved for TVk, and = 
6 k/ fe, to give, 

A^k = (iVb - 1) ^ ; &K = ^ (29) 

Since lu = = ^b/^k> ^nd |ei = substi¬ 

tuting from Eq. (1^ for Nk and 6 ^ into Eqs ([5]) and ([3]) , 
and using Eq. ( 0 ), leads to the following equations for 
the blob scaling variables in terms of bead-spring chain 
parameters. 


where. 


X = (iVb - 1) 


Cl 


36 

6-b5 


Ci = 


9 (TVb - 1)2 6^ 


[N^q^ l*s (6 + 5)2j 


1/3 


(30) 

(31) 


(32) 


It is clear from the expressions above that various com¬ 
binations of bead-spring chain parameters can result in 
identical values for X and Y. Table Hill displays the var¬ 
ious sets of values of bead-spring chain parameters, and 
the resultant sets of scaling variables X and Y, that have 
been used in the simulations reported here, and the cor¬ 
responding bare-model parameter values. Two to three 
different values of iVb are typically used for each value 
of X. Note that Y can be varied independently of X 
by varying (and keeping all other bead-spring chain 
parameters fixed). As a result, the dependence of prop¬ 
erties of dilute polyelectrolyte solutions on X and Y can 
be explored by carrying out BD simulations, and their 
extent of universality assessed. 

In terms of bare-model parameters, the Manning pa¬ 
rameter, 7 o, which is defined as the number of elementary 
charges per Bjerrum length^iS is given by, 70 = fls- The 
OSFKK scaling theory assumes that 70 1, since coun¬ 

terion condensation, which occurs for 70 > 1 , is not taken 
into account— The last column in Table Hill shows that 
the Manning parameter is significantly smaller than one 
in all the simulations reported here. 


D. Mapping from experiments to blob variables 

While the results of the simulations reported here have 
not been compared directly with experimental observa¬ 
tions, it is worth outlining how experimental data on di¬ 
lute polyelectrolyte solutions can be mapped onto blob 
variables, to facilitate future comparisons. In this con¬ 
text, the various equations for evaluating the blob vari¬ 
ables X and Y for an arbitrary polyelectrolyte solution 
are derived in Appendix [B] 

In the current simulation, the behaviour of various 
static properties and the diffusivity of dilute polyelec¬ 
trolyte solutions has been explored for a range of blob 
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(a) 



(b) 


FIG. 2. (Color online) Measures of mean chain size versus the 
reduced screening length, V, at various values of the number 
of blobs, X: (a) scaled end-to-end vector, (b) scaled 

radius of gyration, Simulation data acquired at dif¬ 

ferent values of chain length, Nh, are displayed at each value 
of X, to demonstrate parameter-free data collapse. The vari¬ 
ous combinations of bead-spring chain parameters that corre¬ 
spond to the displayed values of Nh, are shown in Table im 
Error bars are of the order of symbol size. 


variable values, 20 < X < 80 and 0.01 < Y < 100. 
In Appendix E the specific example of an aqueous so¬ 
lution of sodium poly(styrene sulfonate) is considered, 
and it is shown how an appropriate choice of tempera¬ 
ture, polymer molecular weight, degree of ionization per 
chain, and salt and monomer concentration can lead to 
blob variables that span this range. Further, parameter 
values are chosen to ensure that the Manning parameter 
7 o ~ 1 in order to avoid counterion condensation. 


III. RESULTS 

A. Mean chain size as a function of X and Y 

From the OSFKK scaling picture depicted schemat¬ 
ically in Fig. [1] we expect the chain size to be inde¬ 


FIG. 3. (Golor online) Demonstration of the existence of log¬ 
arithmic corrections to scaling at F = 30 (in the blob-pole 
regime) for the scaled end-to-end vector, (a) The dotted line 
shows that departs from the X^ scaling predicted by 

OSFKK theory for X > 30. (b) A power law ht through the 
data for re/(CeiX) versus InX (solid line) shows that the loga¬ 
rithmic corrections have the expected exponent of 1/3 (within 
error bars). The various combinations of bare parameters cor¬ 
responding to the value of Nh are shown in Table imi Error 
bars are of the order of symbol size. 

pendent of Y in the limits of F —>■ 0 (regime VI) and 
F —>■ 00 (regime II) for X ^ 1, while at intermediate 
values of F (in the crossover regimes III, IV, and V), 
we expect the chain size to depend on F. This expecta¬ 
tion has been verified by several prior simulations, which 
have examined the scaling of the end-to-end vector 
Figs. [21(a) and|2](b) show that in the current simulations 
as well, both the scaled end-to-end vector, {r1)/Q^, and 
the scaled radius of gyration, {r'^)/Q\, asymptotically ap¬ 
proach constant values at the two extremes of F. Note 
that the asterisk in has been dropped for simplicity of 
notation, since the non-dimensional character of the blob 
size is clear from the context. In the crossover regime, the 
chain size increases monotonically with F as the chain 
starts to swell due to increasing electrostatic repulsion 
with increasing Debye screening length. As mentioned 
previously, our BD simulations are unable to distinguish 
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(a) 



(b) 


FIG. 4. (Color online) Demonstration of the existence of weak 
logarithmic corrections to scaling at K = 30 (in the blob-pole 
regime) for the scaled radius of gyration, (a) The dotted line 
shows the onset of departure from scaling for at 

large values of X. (b) A power law fit through the last four 
data points for rg/{^eiX) versus InX (solid line) shows that 
the logarithmic corrections have an exponent smaller than 
1/3. The various combinations of bare parameters corre¬ 
sponding to the value of Nh are shown in Table imi Error 
bars are of the order of symbol size. 


between the different regimes in the crossover region due 
to the computational cost of simulating the extremely 
long chains that would be required for such a purpose. 

Notably, as long as the value of X is the same, the 
predicted values of both (Cg)/^gi, and {r^)/^e\^ inde¬ 
pendent of the specific choice of bead-spring chain pa¬ 
rameters, at all values of Y. For instance at X = 48, 
both Nb = 52 and iVb = 60, lead to identical results. 
Such a parameter-free data collapse was previously shown 
to occur for BD simulation predictions of {r1)/Q^ by 
Pattanayek and Prakash^iiwho also noted that the only 
constraint to the simulations appeared to be to ensure 
that there were more beads than blobs in a chain, i.e., 
Nb>X. 

The OSFKK scaling picture ignores the presence of log¬ 
arithmic corrections in the blob-pole regime. However, as 
seen in Appendix and as shown previously by several 


different analytical studiesfid^ such corrections are ex¬ 
pected to arise in this regime, and previous molecular 
simulations have shown that such is indeed the case^i In 
the absence of logarithmic corrections, the scaled end- 
to-end vector, (re)/^gi, is expected to scale as X'^ in 
regime II (see Table IJ). It is clear from Fig. [3] (a) that 
while this scaling is obeyed for values of Ai < 30, there 
is a departure at larger values of AT. A similar devia¬ 
tion from the expected scaling in regime II was observed 
previously in the BD simulations reported by Pattanayek 
and Prakashiii However, they did not verify if the source 
of the deviation was due to logarithmic corrections. In 
terms of the blob scaling variable X, Eq. [5] shows that 
the logarithmic correction term has an exponent of 1/3. 
We expect therefore that a plot of re/(^eiA') versus InAi 
would be a straight line with a slope of 1/3. Fig. [3] (b) 
shows that the logarithmic corrections to OSFKK scal¬ 
ing, picked up by the BD simulations, do in fact have the 
expected dependence on X. 

While the radius of gyration is also a measure of mean 
chain size, Figs. H] reveal that the logarithmic corrections 
to scaling for (rg)/^gi in the blob-pole regime are much 
weaker than for {r"l)/^^y This could be attributed to the 
fact that the radius of gyration not only depends on the 
stretched axis of the chain (which is expected to be in the 
direction of the end-to-end vector), but also depends on 
the dimensions of the chain perpendicular to the elon¬ 
gation axis. The scaling of the chain in these lateral 
dimensions is generally considered to obey random walk 
statistics. We examine the validity of this expectation 
below, when we consider the eigenvalues of the gyration 
tensor. 

In the limit of long chains, the ratio of the mean square 
end-to-end vector to the mean square radius of gyra¬ 
tion, (rg)/(Cg), has a universal value of 6 for ideal chains, 
while it is equal to 12 for rigid rods.— We anticipate, 
therefore, that for the polyelectrolyte solutions consid¬ 
ered here, the value of the ratio will increase from being 
close to 6 in the limit K ^ I (where the chain behaves 
ideally), towards a value of 12, as K —>■ oo, where the 
polyelectrolyte chain adopts a stretched rodlike config¬ 
uration. Fig. [S] (a) shows that the ratio does approach 
a constant value close to 6 as K approaches zero, and 
increases monotonically with increasing Y. However, 
though the ratio levels off as the chain enters the blob- 
pole regime in the limit of large Y, the asymptotic value 
is less than 12, for the values of X considered here. The 
departure from the limiting value for rigid rods can be 
considered to reflect the degree of flexibility remaining in 
the chain in this regime. Interestingly, the ratio appears 
to be independent of X for a range of values of K, start¬ 
ing with the ideal regime, and progressing well into the 
crossover regime (i.e., the red and blue symbols overlap 
in Fig. [3] (a)). The ratio starts to become X dependent 
for values of Y approaching the blob-pole regime. As was 
seen in Figs. [3] (b) and|l](b) above, though both (r^) and 
(Cg) depart from the OSFKK scaling theory predictions 
in this regime due to the presence of logarithmic correc- 
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X 

(b) 

FIG. 5. (Color online) The variation of the ratio {r1) 
with: (a) the reduced screening length V, for two different 
values of the number of blobs X, and (b) the number of elec¬ 
trostatic blobs X, at T = 30. Error bars are of the order of 
symbol size in (a), and are explicitly displayed in (b). 


tions, the strength of these corrections is different in the 
two cases. This difference is responsible for a persistent 
dependence of (?'g)/(^) on X in the blob-pole regime, as 
can be seen in Fig. (b), where the value of the ratio 
is plotted as a function of X, at F = 30. We can an¬ 
ticipate, however, that with increasing values of X (i.e., 
chain length), the appearance of logarithmic corrections, 
and the consequent onset of non-universal behaviour, is 
postponed to larger and larger values of Y, due to a delay 
in the inception of the blob-pole regime (see Fig. [T|). 


B. Eigenvalues of the radius of gyration tensor 

The dependence of the three scaled eigenvalues of the 
radius of gyration tensor, {Xl)/Qy, and (A|)/Cei, 

on the reduced screening length F, is shown in Figs. HI 
As was observed for and (rg)all the eigenval¬ 

ues increase in the crossover regime from constant values 
at small values of F, to constant values in the blob-pole 
regime at large values of F. The collapse of data for dif- 




0.01 0.1 1 10 100 


Y 


(b) 



FIG. 6. (Color online) Scaled eigenvalues of the gyration ten¬ 
sor as a function of F, at two values of X: (a) the smallest 
eigenvalue (Af)/^ei, (b) the intermediate eigenvalue (Ai)/^ei, 
and (c) the largest eigenvalue (Aa)/^^). The insets in (a) 
and (b) display the dependence of the unsealed eigenvalues 
on the non-dimensional Debye length close to the over¬ 
shoot, for two different chain lengths, Nh, at the same value 
of X. Data obtained with the multi-chain algorithm, with 
c/c* = 10~®. Error bars are of the order of symbol size. 


ferent bead-spring chain model parameters, when repre¬ 
sented in terms of blob scaling variables, is also observed 
for these properties. Clearly, the chain size increases with 
increasing electrostatic repulsion between blobs, not only 
in the direction of maximum chain stretching ((A§)), but 
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(a) 



(b) 


FIG. 7. (Color online) Demonstration of the existence of log¬ 
arithmic corrections to scaling at y = 30 (in the blob-pole 
regime) for the largest scaled eigenvalue of the gyration ten¬ 
sor. (a) The dotted line shows that {^)/^^\ departs from 
the scaling predicted by OSFKK theory for X > 30. (b) 
A power law fit through the the last four data points for 
{A|)/^eiA versus InX (solid line) shows that the logarith¬ 
mic corrections have an exponent close to the expected value 
of 1/3. Data was obtained using the multi-chain algorithm 
with c/c* = 10“®. Error bars are of the order of symbol size 
in (a), and are explicitly displayed in (b). 


also in the lateral directions ((Af) and (A 2 )). The sig¬ 
nificant differences in the magnitude of the chain dimen¬ 
sions in the three different principal directions clearly 
illustrates the highly anisotropic shape of the chain in all 
regions of the phase diagram. This is examined in more 
detail in terms of shape functions in the section below. A 
curious observation, for which we do not have an obvious 
explanation, is the presence of an overshoot in (Ai)/5g[ 
and {^ 2 )/Q\ foi' values of Y at the threshold of the blob- 
pole regime, i.e. just before the chain enters the final 
screening-length-independent regime from the crossover 
regime. The continued stretching of (A|) with increas¬ 
ing Y at the verge of the blob-pole regime appears to 
be accommodated by a shrinkage in (Af) and (A|). The 
overshoot also becomes more pronounced with increas¬ 
ing X. The insets in Figs. [5] (a) and (b) demonstrate 
that even the occurrence of an overshoot, which is dis- 



X 


(a) 



FIG. 8. (Golor online) Dependence on the number of blobs 
X, of (a) the smallest scaled eigenvalue of the gyration tensor, 
(Ai)/Cei) and (b) the intermediate scaled eigenvalue of the 
gyration tensor, (A^j/CeU at a reduced screening length, Y = 
30. Data was obtained using the multi-chain algorithm with 
c/c* = 10“®. Error bars are of the order of symbol size. 


tinct for different bead-spring chain parameters, collapses 
onto a single curve when represented in terms of scaling 
variables—signifying that blob variables accurately cap¬ 
ture the essential physics, even for phenomena that go 
beyond the simple picture upon which they were origi¬ 
nally based. The use of finitely extensible springs in the 
simulations is probably not responsible for the overshoot, 
and it is more likely to be an electrostatic phenomenon, 
since the springs are weakly stretched even in the blob- 
pole regime. For instance, from the value of {r1)/Q^ in 
Fig. m (a) ((r 2)/^2 _ 2000, for = 60 and X = 48, 
in the limit of large T), and the values of ^*1 and L* 
corresponding to this value of N\y and X in Table IIIII 
(^*j = 1.85 and L* = 471), we can see that the ratio of 
the end-to-end vector to the contour length is less than 
16%. 

The scaling of with X in the blob-pole regime was 
seen in Fig.[3](b) to exhibit logarithmic corrections in line 
with the prediction of Eq. ®. Since the largest eigen¬ 
value of the radius of gyration tensor is expected to be 
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TABLE IV. Normalized eigenvalues of the radius of gyration tensor for neutral random walk polymers compared with those 
for a polyelectrolyte chain (PE) in the ideal chain regime of the phase diagram. Predictions by Koyama^^ and Wei— were 
obtained from approximate analytical calculations for infinite Gaussian chains, while those by Sole—, Sole and Stockmayer— 
and Kranbuehl and Verdier— were obtained with lattice Monte Carlo simulations. Values reported from Theodorou and 
Suter— are for atactic polypropylene, obtained from Monte Carlo simulations based on a rotational isomeric state model, 
while those from Zifferer— are from lattice Monte Carlo simulations of “nonreversal random walks”, which are considered to 
represent 6-solvents. Data for the polyelectrolyte solution was obtained using the multi-chain algorithm, with c/c* = 10~®. 


Solution 

Chain length 

(A^)/(r■^) 



(A^)/(Af> 

(Ai>/(A() 

Neutrali^ 

Ab —>■ oo 

0.0646 

0.175 

0.754 

2.5 

10.6 

Neutrali^ 

Ab —>■ oo 

0.071 

0.179 

0.750 


10.554 

hJeutraSSS 

Ab = 100 

0.065 

0.175 

0.75 

2.7 

11.7 

NeutraU^ 

Ab = 63 

0.065 

0.176 

0.76 

2.70 

11.7 

Neutrali^ 

Ab = 999 

0.060 ± 0.004 

0.17 ±0.01 

0.77 ±0.06 

2.9 ±0.2 

12.7 ±0.9 

Neutral^ 

Ab —>■ oo 

0.0633 

0.1721 

0.7645 



PE {Y = 0.01) 

to 

0.0631 ± 0.0003 

0.1734 ±0.0009 

0.7626 ± 0.0061 

2.75 ±0.02 

12.10 ±0.12 


TABLE V. Normalized eigenvalues of the radius of gyration tensor for various values of the number of blobs X, at a reduced 
screening length, Y = 30. Data was obtained using the multi-chain algorithm, with c/c* = 10“®. 


A 

{ K )/ K ) 



(A|)/(Af> 


1 - (A0/(A|> 

20.3 

0.0168 ±0.0001 

0.0734 ± 0.0003 

0.9099 ± 0.0026 

4.37 ±0.02 

54.2 ±0.22 

0.771 ±0.001 

24 

0.0148 ± 0.0001 

0.0666 ± 0.0004 

0.9186 ± 0.0038 

4.52 ±0.03 

62.3 ± 0.30 

0.779 ± 0.002 

27.6 

0.0136 ±0.0001 

0.0636 ±0.0011 

0.9228 ± 0.0078 

4.69 ±0.09 

68.1 ±0.88 

0.787 ± 0.004 

30.7 

0.0125 ±0.0001 

0.0578 ± 0.0004 

0.9297 ± 0.0034 

4.63 ±0.04 

74.5 ± 0.43 

0.784 ± 0.002 

34.8 

0.0114 ±0.0001 

0.0545 ±0.0011 

0.9341 ± 0.0076 

4.77 ±0.11 

81.8 ± 1.12 

0.790 ± 0.004 

37.7 

0.0108 ±0.0001 

0.0536 ±0.0004 

0.9356 ± 0.0030 

4.95 ±0.04 

86.5 ± 0.40 

0.798 ± 0.001 

48 

0.0094 ± 0.0002 

0.0488 ±0.0013 

0.9418 ± 0.0074 

5.22 ±0.16 

100.7 ± 1.7 

0.808 ± 0.006 


in the direction of maximum chain stretching, we expect 
\/ (A|)/^ei to also exhibit logarithmic corrections accord¬ 
ing to Eq. ([5|). Figure [7] (a) displays the departure from 
the OSFKK scaling prediction in the blob-pole regime, 
while Fig. [7] (b) shows that the logarithmic corrections 
have an exponent close to the expected value of 1/3. 

According to OSFKK theory, the scaling of chain di¬ 
mensions lateral to the elongation axis is expected to re¬ 
main unperturbed in the blob-pole regime, i.e., to obey 
random walk statistics. To our knowledge, there are no 
theories that describe modifications to this scaling be¬ 
haviour due to the nonuniform stretching of the chain 
(unlike the refined theories developed for scaling in the 
direction of maximum stretching). The dependence of 
chain dimensions on the number of blobs, in the direc¬ 
tions perpendicular to the stretching direction, is exam¬ 
ined in Figs. [5] (a) and (b), where, (Af)/^gj and (A|)/5g[ 
are plotted as functions of X. For the limited range of 
values of X examined here, the smallest and intermedi¬ 
ate eigenvalues appear to obey power law scaling with 
exponents that lie somewhere between ideal chain and 
blob-pole scaling. 

Eigenvalues of the radius of gyration tensor for neutral 
polymer chains are usually reported in terms of ratios, 
either between individual eigenvalues, or with the mean 
square radius of gyration, since they are expected to at¬ 
tain universal values in the limit of long chains. For a 


chain with a spherically symmetric shape about the cen¬ 
tre of mass, we expect (A^)/(rg) = 1/3, for i = 1,2,3, 
and (Af)/(A|) = 1 for all combinations i and j. Pre¬ 
dictions reported previously in the literature for flexible 
neutral chains in d-solutions, from a variety of differ¬ 
ent approaches, are displayed in Table IIVI The strong 
asymmetry in the shapes of chains is clearly apparent 
from the extent of departure from the expected values. 
Steinhauser— has carried out molecular dynamics sim¬ 
ulations to examine the influence of solvent quality on 
the ratios of eigenvalues, and noted that while their val¬ 
ues are relatively similar in good and theta solvents, 
they change significantly in poor solvents, with values 
that imply that chains become much more spherical as 
they collapse into globules. From OSFKK theory, we ex¬ 
pect polyelectrolyte chains to behave like ideal chains for 
K <C 1. As can be seen from the last row of Table El 
this is indeed the case for the results obtained with the 
current simulations, with all the different ratios being in 
good agreement with results for neutral chains. 

The dependence of the eigenvalues scaled with the 
mean square radius of gyration, on the reduced screen¬ 
ing length Y, is displayed in Figs. [9] Starting with the 
limiting values reported in Table IIVI for polyelectrolyte 
chains in the ideal chain regime, both (Af)/(rg) and 
(not shown here) decrease with increasing Y 
in the crossover regime until they reach their asymp- 





































13 



Y 



0.01 0.1 1 10 100 
V 

(b) 

FIG. 9. (Color online) Dependence of the eigenvalues of the 
radius of gyration tensor, normalised by the radius of gyra¬ 
tion, on the reduced screening length Y, at two values of the 
number of blobs, X: (a) the smallest eigenvalue (Ai)/(rg), 
and (b) the largest eigenvalue {^)/{r^)- Data was obtained 
using the multi-chain algorithm, with c/c* = 10“®. Error 
bars are of the order of symbol size. 


totic values in the blob-pole regime. The reason for 
the decrease in the magnitude of these ratios is because 
(j'g) = (Af) -b (A|) -f (Ag) is dominated by the behaviour 
of (A§), which can be seen from Figs. [6] to increase much 
more rapidly with Y than both (Af) and (A^). For the 
same reason, the ratio (A§)/(rg) increases with Y, nearly 
approaching a value of 1 in the limit of large Y. The over¬ 
shoot in (Ag) and (A|) that was observed in Figs. [6] (a) 
and (b) at the threshold of the blob-pole regime, is not 
noticeable when the eigenvalues are scaled with (r|) in 
place of (as displayed in Fig. |n](a) for {X\)/{r‘^)). 

The extent of asymmetry in chain shape in the plane 
perpendicular to the elongation axis can be gauged by 
the quantity, 1 — (Af)/(A 2 ), which would be zero if the 
shape was symmetric. For neutral chains, and for poly¬ 
electrolyte chains in the ideal chain regime, using the 
values of (Af)/(rg) and (A|)/(rg) reported in Table HVl 
we can calculate 1 —(Af)/(A|) ~ 0.63. Fig.fTUlshows that 
with increasing electrostatic repulsion between blobs, the 
already highly asymmetric shape becomes even more so, 



FIG. 10. (Color online) Asymmetry in chain shape in the 
plane perpendicular to the stretching axis as a function of the 
reduced screening length, Y. Data was obtained using the 
multi-chain algorithm, with c/c* = 10~®. 


levelling off to a constant value in the limit of large Y. 

A notable aspect of the normalised eigenvalue ratios 
displayed in Figs. [9] and [10] is that their values are inde¬ 
pendent of the number of blobs X in the chain, for a wide 
variety of values of Y, ranging from the ideal chain regime 
and well into the crossover regime. The ratios become 
X dependent only as the chains approach the blob-pole 
regime with increasing screening length. This suggests 
that for sufficiently long chains, as in the case of neutral 
polymer chains, the normalised eigenvalues of polyelec¬ 
trolyte chains attain universal values in those parts of 
the phase diagram that he outside the blob-pole regime. 
Within the blob-pole regime, we expect that ratios of 
quantities that scale nearly identically with X will attain 
roughly constant values, such as (A|)/(rg) and (A|)/(rg) 
(which follows from the dominance of the growth of (A|) 
with X, compared to (Af) and (A|)), while ratios of quan¬ 
tities that do not scale identically, would retain a de¬ 
pendence on the number of blobs. This is examined in 
Table El where the dependence of a variety of ratios on 
the number of blobs X, at the particular value, Y = 30, 
which is in the blob-pole regime, is tabulated. At any 
fixed value of Y, however, we expect chains to move out 
of the blob-pole regime as X —>■ oo, and for all the ratios 
to become universal. 


C. Shape functions 

For chain shapes with tetrahedral or greater symme¬ 
try, the asphericity B = 0, otherwise B > 0. For 
chain shapes with cylindrical symmetry, the acylindric- 
ity C = 0, otherwise C > 0. With regard to the degree 
of prolateness, its sign determines whether chain shapes 
are preponderantly oblate {S,S* G [—0.25,0]) or prolate 
{S,S* G [0,2]). The relative anisotropy (k^ and k^*), 
on the other hand, lies between 0 (for spheres) and 1 
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(c) 


(d) 


FIG. 11. (Color online) The dependence of shape functions on the reduced screen length Y, at two values of X: (a) the 
normalised asphericity B/{rg) (see Eq. (1191) '). (b) the normalised acylindricity C/ (rg) (see Eq. (1201) 1. (c) the degree of prolateness, 
S* (see Eq. (1211) 1. and, (d) the shape anisotropy, k^* (see Eq. (1231) 1. When not displayed explicitly, error bars are of the order 
of symbol size. Data obtained with the multi-chain code, with c/c* = 10~®. 


TABLE VI. Shape functions for neutral random walk polymers compared with those for a polyelectrolyte chain (PE) in the 
ideal chain regime of the phase diagram. Predictions by Steinhauser— were obtained with molecular dynamics simulations, 
while those by Bishop and Michels— and Bishop et al.— were obtained with Brownian dynamics, and off-lattice Monte Carlo 
simulations, respectively. Methods used to evaluate shape functions in Refs, [2^ and[23 are given in the caption to Table ITVl 
Data for the polyelectrolyte solution was obtained using the multi-chain algorithm, with c/c* = 10~®. 


Solution 

Chain length 

B/K) 

c/K) 

s* 

S 



Neutral— 

Nh oo 



0.4750 

0.8875 

0.3943 

0.5263 

NeutralSi 

Nb = 999 

0.66 ± 0.04 

0.106 ± 0.007 



0.41 ± 0.01 


Neutral^ 

OO 

II 



0.460 ± 0.021 

0.859 ± 0.245 

0.384 ± 0.011 

0.518 ± 0.095 

Neutrals 

Nb -1 oo 





0.397 ± 0.001 

0.529 ± 0.001 

Neutral;^ 

Nb —>• oo 

0.6255 ± 0.0005 




0.394 ± 0.003 

0.5221 ± 0.0013 

Neutrals 

Nb -t oo 



0.4751 

0.8875 

0.3943 

0.5264 

PE (Y = 0.01) 

II 

to 

0.6453 ± 0.0056 

0.1103 ± 0.0006 

0.4751 ± 0.0024 

0.8810 ± 0.0205 

0.3955 ± 0.0011 

0.5253 ± 0.0074 


(for rods). In this context, the values of the normalised 
shape functions B/(r^) and C/{r'^), S and S*, and and 
*, predicted by various techniques for neutral polymer 
chains are compared with the predictions of the current 
simulations for polyelectrolyte chains in the ideal chain 
regime, in Table IVII Clearly, in this regime, polyelec¬ 
trolyte chains behave identically to neutral random walk 


chains, and share the same high degree of anisotropy. 

The dependence of the shape functions B/{r‘^), C/ (r^), 
S* and on T, for two different values of X, is dis¬ 
played in Figs. [TT] As we might expect, chains appear to 
become more aspherical, more cylindrical, more prolate, 
and more rodlike, with increasing electrostatic repulsion 
between the blobs. The functions S and behave simi- 
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(a) 




(b) 



(c) 


(d) 


FIG. 12. (Color online) The dependence of shape functions on the number of blobs X, in the blob-pole regime, at a particular 
value of the reduced screening length Y = 30: (a) the scaled asphericity (b) the scaled acylindricity C/^eii (c) the degree 

of prolateness S, and, (d) the shape anisotropy Error bars are of the order of symbol size in (a) and (b). Data obtained 
with the multi-chain code, with c/c* = 10“®. 


larly to S* and k^*. 

The independence of the values of the shape functions 
from the number of blobs X, in all the regimes except 
the blob-pole regime, suggests that for sufficiently long 
chains, at any given value of the reduced screening length 
Y, the shapes of polyelectrolyte chains are universal. In 
the blob-pole regime however, as noted earlier, the ap¬ 
pearance of corrections to scaling (to varying extents in 
the three principal directions), leads to a dependence of 
chain shapes on the chain length. The nature of this 
dependence on X, at a particular value of T = 30, is 
displayed in Figs. [T^] for i?/Ceu S and The 

functions S* and behave similarly to S and . 


D. Translational diflusivity interpreted with the OSFKK 
scaling picture 


they are functions of Y in the crossover regime between 
these two limits. By defining the Zimin diffusivity of an 
electrostatic blob through the expression, 

" 67ri7,(eei/2) 

we hnd that a similar pattern is obeyed by the ratio 
D/D^, where D is the translational diffusivity of the 
polyelectrolyte chain as a whole. 

In the ideal chain regime, since D ^ k^T/(dTrrjsre), we 
expect D/D^ ~ In the blob-pole regime, an ex¬ 

pression for the diffusivity can be derived by drawing an 
analogy with the shish-kebab model for rodlike polymers. 
In the latter, a rodlike polymer of length L is modelled as 
N = (L/d) beads of diameter d, placed along a straight 
line. The translational diffusivity of such a shish-kebab 
can be shown to be;^ 


The dependence on the reduced screening length Y, 
of all the static properties examined here so far, has fol¬ 
lowed a similar pattern: they exhibit Y independent be¬ 
haviour in the ideal chain and blob-pole regimes, while 


D = 


iTTrjsL 


(34) 


By mapping d —>■ ^ei, and L —^ re, it follows that the 
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FIG. 13. (Color online) Dependence of the ratio of the trans¬ 
lational diffusivity of the chain D, to the diffusivity of the 
electrostatic blob D^, on the reduced screening length Y, at 
two values of the number of blobs, X. Error bars are of the 
order of symbol size. 


diffusivity of a blob-pole is given by, 


D = 


In(^) 

37r?7sre 


ksT 


(35) 


Within the OSFKK scaling ansatz, since (rg/^ei) ~ X in 
the blob-pole regime, we see that. 


D _ ln(X) 


(36) 


If logarithmic corrections to the scaling of re (according 
to Eq. ([5l)) are taken into account, we expect, 


D [ln(X)]^/^ 
~ X 


In either case, we see that the ratio D/D^ is independent 
of y in the both the limiting regimes of small and large Y. 
In between these two limits, we expect D/D^ to depend 
on Y. 

These arguments are substantiated in the plot of D/D^ 
versus T, at two values of X, displayed in Fig. [131 We see 
that for each value of X, the ratio D/D^ crosses over from 
a constant value in the ideal chain regime to a constant 
value in the blob-pole regime. As has been observed in 
the case of all the properties examined so far, the use 
of blob scaling variables to interpret the results of BD 
simulations leads to behaviour that is independent of the 
specific choice of bead-spring chain parameters. 

The importance of accounting for logarithmic correc¬ 
tions to the scaling of with X in Eq. (1551) . is stud¬ 
ied in Figs, [m The plot of D/(D^lnX) versus 1/X 
in Fig. [14] (a) shows that while the scaling described by 
Eq. ((5S1) is obeyed at relatively small values of X, there 
is a departure from OSFKK scaling for X > 30. On the 
other hand, the plot oi D/{D^X~^) as a function of InX 


a 




(a) 



FIG. 14. (Color online) Demonstration of the importance 
of accounting for logarithmic corrections to the scaling of Ve 
with X in Eq. (1351) for the translational diffusivity D in the 
blob-pole regime: (a) Departure from OSFKK scaling at large 
values of X (see Eq. (1361) 1 (b) Scaling behaviour in agreement 
with Eq. (EZl), with the expected exponent of InX « 2/3. 
Error bars are of the order of symbol size. 


displayed in Fig. [T3] (b) shows that Eq. (1571) is indeed 
obeyed, with the exponent of In X close to the expected 
value of 2/3. 

A comparison of the magnitude of the diffusivity in the 
ideal chain and blob-pole regimes in Fig.[T31 suggests that 
strongly screened chains {Y —)> 0) diffuse more than twice 
as fast as nearly unscreened chains {Y :s> 1), and that 
this ratio increases with X. This behaviour is examined 
in more detail in Fig. [T51 which displays the ratio of the 
diffusivity of a chain with Y = 0.01 to that of a chain 
with y = 30, for various values of X. When logarithmic 
scaling corrections are ignored, we expect this ratio to 
scale as X 2 . As can be seen from Fig. [151 this scaling is 
approximately satisfied. As in other examples considered 
here, the inclusion of logarithmic corrections in the blob- 
pole regime leads to a better fit of the data, as shown in 
the inset to Fig. [T51 

From the scaling of re/^ei (see Eq. ([S])) and D/D^ (see 
Eq. (1571) ) with X in the blob-pole regime, it can be seen 
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FIG. 15. (Color online) Dependence on the number of blobs 
X, of the ratio of the translational diffusivity of the chain in 
the ideal chain regime {Y = 0.01) to the diffusivity of a chain 
in the blob-pole regime {Y = 30). Inset shows the importance 
of accounting for logarithmic corrections to scaling. Error 
bars are of the order of symbol size. 


that the ratio, 


{rl)/a ^ {rl)/D ^ ^3 

D/D(^ Tj 

where, r^, the blob relaxation time, is defined by, rj = 
Notably, this ratio is free from logarithmic cor¬ 
rections even in the blob-pole regime. Figure [16] (a) dis¬ 
plays the dependence on the number of blobs X, of the 
ratio of the relaxation time in the blob-pole regime to the 
relaxation time of an electrostatic blob, t/tj. Though 
the relaxation time r has been dehned in terms of {r^) 
instead of (r^) (see Eq. (l26)l ). we see from Fig.|T6|(a) that 
(in the range of X values that have been considered), any 
residual logarithmic corrections are extremely weak, and 
that the ratio scales with X with the expected power law 
exponent. 

The difference in the conformations of the chain in the 
two regimes, leads to the long-time relaxation in the blob- 
pole regime being nearly an order of magnitude slower 
than in the ideal chain regime, as displayed in Fig.fTOKb). 
As can be seen from the figure, the ratio of the two lim¬ 
iting relaxation times, 'r(v= 30 )/'''(i'=o.oi)) increases with 
the number of blobs X with the expected power law, X i. 

The ratio C/rd, defined by. 

C/rd = ^ (38) 

where, rh is the hydrodynamic radius, rn = 
U'bT/(Q'KT jsD), has an universal value for neutral poly¬ 
mer solutions in the long chain limit, since both Tg and 
Th scale identically with chain length i^"^d^'^° In partic¬ 
ular, by extrapolating finite chain data acquired from 
highly accurate BD simulations, to the long chain limit, 
Sunthar and Prakash"^"*^ have shown that for 0-solutions, 



(a) 



X 


(b) 


FIG. 16. (Color online) Dependence on the number of blobs 
X, of (a) the ratio of the relaxation time of the chain in the 
blob-pole regime to the relaxation time of an electrostatic 
blob, and (b) the ratio of the relaxation time of the chain in 
the blob-pole regime (at Y = 30) to the relaxation time in 
the ideal chain regime (at Y = 0.01). Error bars are of the 
order of symbol size. 

C/^D = 1-38 ± 0.01. We anticipate that for the polyelec¬ 
trolyte solutions considered here. C/rd will have a similar 
value in the ideal chain regime. 

In the blob-pole regime, by substituting for D from 
Eq. (l35l) into Eq. (l38)) . we see that, 

C/rdI n A (39) 

As a result. C/rd does not have a universal value in this 
regime, but rather depends on the number of blobs X. 
In the context of the current simulations, Figs. |3| (b) 
and|4|(b) show that the scaling expressions, re/(5eiA) ~ 
(lnA)°-^®, and rg/(^eiA) ^ (lnA)°-^^, respectively, pro¬ 
vide good fits to the simulation data. As a result, we 
expect from Eq. (l3^ that. C/rd (lnA)°-®®. 

EigureflTl (a) displays the dependence of C/rd on the 
reduced screening length F, for two different values of 
X. We see that for small values of Y, the ratio has a 
constant universal value which is close to that for neutral 
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(a) 



InX 

(b) 

FIG. 17. (Color online) Ratio of the radius of gyration to the 
hydrodynamic radius, I/rd (see Eq. ([38])), as a function of 
(a) the reduced screening length Y, at two different values of 
the number of blobs X, and (b) demonstration of logarithmic 
corrections to scaling in the blob-pole regime, at a constant 
value of the reduced screening length, Y = 30. Error bars are 
of the order of symbol size. 


chains. The ratio increases, while remaining universal 
(i.e., independent of X) for increasing values of Y in the 
crossover regime, before levelling off to a non-universal 
value in the blob-pole regime. The dependence of C/rd 
on InX in the blob-pole regime, at the particular value 
Y = 30, is displayed in Fig. |T7] (b). We see that the 
exponent of In X is indeed close to the expected value of 
0 . 86 . 


IV. CONCLUSIONS 

The size, shape and diffusivity of a weakly-charged 
polyelectrolyte chain in solution have been examined in 
the limit of low polymer concentration using Brown¬ 
ian dynamics simulations of a coarse-grained bead-spring 
chain model, with Debye-Hiickel electrostatic interac¬ 
tions between the beads. Simulation results have been 
recast in terms of the scaling variables X (the num¬ 


ber of electrostatic blobs), and Y (the reduced screening 
length), which are defined within the framework of the 
OSFKK blob scaling theory. While the root-mean-square 
end-to-end vector and radius of gyration have been used 
to represent the mean size of the chain, various func¬ 
tions defined in terms of the eigenvalues of the radius 
of gyration tensor have been used to characterise chain 
shape in all parts of the {X,Y} phase-space. The trans¬ 
lational diffusion coefficient of the chain, which is a dy¬ 
namic property, has been determined accurately from the 
displacement of the chain centre of mass, by taking hy¬ 
drodynamic interactions into account through incorpora¬ 
tion of the Rotne-Prager-Yamakawa tensor into the BD 
algorithm. 

The key results of the present work are summarised in 
the list below: 

1. Interpretation of simulation results in terms of blob 
scaling variables X and Y leads to a description 
of solution behaviour independent of the level of 
coarse-graining, i.e., of the choice of the number of 
beads N\^, and other parameters in the bead-spring 
chain model. This should prove extremely useful 
for comparing simulation results with experiments. 
The procedure for mapping from experimental vari¬ 
ables to blob scaling variables, in order to facilitate 
this comparison, has been provided in Appendix IbI 

2. Three broad domains of behaviour can be clearly 
identified: (i) the ideal chain regime corresponding 
to small values of Y, (ii) the crossover regime cor¬ 
responding to intermediate values of Y, and (iii) 
the blob-pole regime corresponding to large values 
of Y. 

3. In the blob-pole regime, BD simulations appear to 
pick up the existence of logarithmic corrections for 
fairly short chains, which are in good agreement 
with predictions of refined scaling theories. 

4. Various universal ratios of eigenvalues, the as- 
phericity, the acylindricity, the degree of prolate¬ 
ness and the shape anisotropy of polyelectrolyte 
chains in the ideal chain regime, have been com¬ 
pared with published results in the literature for 
flexible neutral random walk polymers. 

5. When suitably scaled, the size, shape and diffu- 
sivity of a chain are independent of the reduced 
screening length Y in the ideal chain and blob-pole 
regimes, and depend on Y only in the crossover 
regime. 

6. The translational diffusivity of the chain in the 
blob-pole regime can be described by drawing an 
analogy with the translational diffusivity of a rod¬ 
like polymer modelled as a shish-kebab. The nor¬ 
malisation parameter that enables the collapse of 
BD data for the translational diffusivity is the 
Zimm diffusivity of the electrostatic blob. 
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7. All properties of polyelectrolyte solutions, when 
suitably normalised, appear to exhibit universal be¬ 
haviour, independent of the number of electrostatic 
blobs X in the chain, in all regimes of the {X, Y} 
phase-space, except in the blob-pole regime, where 
the occurrence of logarithmic corrections to scaling 
leads to non-universal behaviour. 
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Appendix A: Scaling of the end-to-end vector in the 
blob-pole regime 


The scaling of the end-to-end vector in the blob-pole 
regime can be derived from a Flory type energy min¬ 
imisation argument. The essential assumption is that 
chains adopt an ellipsoidal shape, with electrostatic inter¬ 
actions causing stretching in the direction of the major- 
axis, while leaving chain conformations in the direction 
of the minor-axes unaffected. As a result, the long half¬ 
axis of the ellipsoid scales as 7’e/2, while the aspect ratio 
scales as The total energy of such a chain 

can be obtained by combining the electrostatic energy of 
an ellipsoid with the elastic energy of a stretched chain. 
Ignoring prefactors of order unity, the electrostatic en¬ 
ergy of an ellipsoid is,— 




MNl 


while the elastic energy of a stretched chain isfi 


fcsT 



From Flory theory, it follows that the equilibrium end-to- 
end vector can be derived by minimising the total energy. 


JL ( \ 

k^T ^ \blNK) 


Ib ifNKf 

Ve 


In 


.bKNi, 


(Al) 


It is convenient to proceed in two steps—first deriving the 
equilibrium end-to-end vector by neglecting the logarith¬ 
mic term in the electrostatic energy, and then accounting 
for it in the next step. Minimising the sum. 


) Te 


with respect to re, leads to, 

Te = Ib bK (A2) 

This is identical to the expression derived by OSFKK the¬ 
ory for the end-to-end vector in the blob-pole regime (see 
Table m. We now assume that including the logarithmic 
term in the electrostatic energy leads to a modification 
of the equilibrium end-to-end vector, 

re = fi l^NKbKX (A3) 


where, the factor x remains to be determined. Substi¬ 
tuting the modified expression for from Eq. (IMI) into 
Eq. (lAll) (and absorbing the common factor A^k(/^ ^b)^ 
into the energy), leads to the following sum. 


- In 
X 


Nl /§ tl X 


which must be minimised with respect to x in order to 
determine x. This leads to. 


=ln 



-I- In a; — 1 


Assuming a; ^ 1, implies. 


ri 

r 4 - 2 -| 


NKpli 


1 . 

3 


(A4) 


(A5) 


where, the factor (1/2) in front of the logarithmic term 
has been introduced for future convenience. The assump¬ 
tion that X is large can be seen to be justified for Vk ^ 1. 

Substituting the expression for x from Eq. into 

Eq. (IA3I) . shows that this two-step procedure leads to a 
prediction of the equilibrium end-to-end vector scaling 
with chain size according to, 

^ ^ [| Vk {in [Nk f|] } ' (A6) 

Clearly, the logarithmic correction to the OSFKK scaling 
expression in the blob-pole regime comes from taking the 
logarithmic term in the expression for the electrostatic 
energy of an ellipsoid into account. Equation (IA6I) is sim¬ 
ilar to the expression derived previously by Dobrynin and 
Rubinstein— This expression takes a particularly sim¬ 
ple form in terms of blob scaling variables (see Eqs. (|1]) 
and (HI)), 

^-A[lnA]3 (A7) 


Appendix B: Mapping experiments to blob variables for 
NaPSS 

Consider a polyelectrolyte chain of molecular weight 
M, dissolved in a solvent at temperature T, with a tem¬ 
perature dependent relative permittivity er(T'). If the 
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TABLE VII. Mapping experimental data for sodium poly(styrene sulfonate) to the blob model 


M = 20600 (Nm = 100, Nk = 10) 


Solution properties 

Blob variables 

«(/) 

cf (I) (mol/ltr) 

T°C 

Zb (nm) 

Zd (nm) 

70 

?el (nm) 

X 

Y 

0.01 (0.1) 

0 (5 X 10“®) 

15 

0.7077 

4321.0 

0.0283 

17.67 

0.2001 

244.5084 

25 

0.7158 

4297.0 

0.0286 

17.61 

0.2016 

244.0463 

35 

0.7247 

4270.0 

0.0290 

17.53 

0.2033 

243.5422 

0.01 (0.01) 

15 

0.7077 

3.056 

0.0283 

17.67 

0.2001 

0.1729 

25 

0.7158 

3.038 

0.0286 

17.61 

0.2016 

0.1726 

35 

0.7247 

3.019 

0.0290 

17.53 

0.2033 

0.1722 

0.5 (0.5) 

15 

0.7077 

0.4.321 

0.0283 

17.67 

0.2001 

0.0245 

25 

0.7158 

0.4297 

0.0286 

17.61 

0.2016 

0.0244 

35 

0.7247 

0.4270 

0.0290 

17.53 

0.2033 

0.0244 

0.3 (3.0) 

0 (1.5 X 10-'^) 

15 

0.7077 

788.9 

0.8492 

1.830 

18.6534 

431.0033 

25 

0.7158 

784.5 

0.8589 

1.824 

18.7951 

430.1889 

35 

0.7247 

779.6 

0.8696 

1.816 

18.9512 

429.3002 

0.01 (0.01) 

15 

0.7077 

3.056 

0.8492 

1.830 

18.6534 

1.6693 

25 

0.7158 

3.038 

0.8589 

1.824 

18.7951 

1.6661 

35 

0.7247 

3.019 

0.8696 

1.816 

18.9512 

1.6627 

0.5 (0.5) 

15 

0.7077 

0.4321 

0.8492 

1.830 

18.6534 

0.2361 

25 

0.7158 

0.4297 

0.8589 

1.824 

18.7951 

0.2356 

35 

0.7247 

0.4270 

0.8696 

1.816 

18.9512 

0.2351 


M = 206000 {Nm = 1000, Nk = 100) 


«(/) 

c™ (I) (mol/ltr) 

T°C 

Zb (nm) 

Zd (nm) 

70 

?el (nm) 

X 

Y 

0.1 (1.0) 

4 X 10“® (4.05 X 10“®) 

15 

0.7077 

151.8 

0.2831 

3.808 

43.1119 

39.8766 

25 

0.7158 

151.0 

0.2863 

3.793 

43.4393 

39.8012 

35 

0.7247 

150.0 

0.2899 

3.777 

43.8001 

39.7190 

0.1 (0.1) 

15 

0.7077 

0.9662 

0.2831 

3.808 

43.1119 

0.2538 

25 

0.7158 

0.9608 

0.2863 

3.793 

43.4393 

0.2533 

35 

0.7247 

0.9548 

0.2899 

3.777 

43.8001 

0.2528 

0.8 (0.8) 

15 

0.7077 

0.3416 

0.2831 

3.808 

43.1119 

0.0897 

25 

0.7158 

0.3397 

0.2863 

3.793 

43.4393 

0.0896 

35 

0.7247 

0.3376 

0.2899 

3.777 

43.8001 

0.0894 

0.35 (3.5) 

4 X 10“® (4.175 X 10-®) 

15 

0.7077 

149.5 

0.9907 

1.652 

229.0977 

90.5376 

25 

0.7158 

148.7 

1.0021 

1.645 

230.8376 

90.3666 

35 

0.7247 

147.8 

1.0146 

1.639 

232.7550 

90.1799 

0.1 (0.1) 

15 

0.7077 

0.9662 

0.9907 

1.652 

229.0977 

0.5850 

25 

0.7158 

0.9608 

1.0021 

1.645 

230.8376 

0.5839 

35 

0.7247 

0.9548 

1.0146 

1.639 

232.7550 

0.5827 

0.8 (0.8) 

15 

0.7077 

0.3416 

0.9907 

1.652 

229.0977 

0.2068 

25 

0.7158 

0.3397 

1.0021 

1.645 

230.8376 

0.2064 

35 

0.7247 

0.3376 

1.0146 

1.639 

232.7550 

0.2060 


monomer concentration in molar units is c™, and the 
degree of ionization per chain is a, then the number of 
counterions in solution per unit volume is, a c“(AiAl000), 
where TVa is Avagadro’s number. If the added salt 
with molar concentration c“, dissociates in solution ac¬ 
cording to, 

AjjjBjj —>■ mA^’’" + nB^~ 

where, z+ and are the cation and anion valences, re¬ 
spectively, then the number of cations A per unit volume 
in solution is c+ = m c™(A^a 1000), and the number of an¬ 
ions B per unit volume is c_ = n c“(iVAl000). It follows 
that the ionic strength I of the solution is given by, 

I = - {rnz\ + nz^) c“(A^a1000) + -z^a c™(IVa1000) 

(Bl) 

where, Zp is the counterion valence. 

The Bjerrum length 1 b can be calculated from Eq. o. 


while the Debye length is given by. 

Id = (B2) 

Mapping the experimental system onto blob variables 
is straightforward once it is mapped onto the bare-model 
parameters {A^k, /, Ib, Id}- In order to do so, it is neces¬ 
sary to know the intrinsic persistence length of the poly¬ 
electrolyte chain, £o, the molecular weight of a monomer, 
Mm, and the monomer length 6m- With this informa¬ 
tion, we can find, (i) the Kuhn step length, 6 k = 2£o, 

(ii) the number of monomers in a chain, = M/Mm, 

(iii) the number of Kuhn steps, IVk = (-/Vm6m)/6K, and 

(iv) the number of charges per Kuhn step, / = aA^m.K, 
where, A^m,K = 6K/6m, is the number of monomers per 
Kuhn step. It is also straight forward to find 1 b = ^b/6k, 
and Id = 1d/6k, at the given temperature and ionic 
strength. Once the bare-model parameters are known, 
the blob scaling variables, {A, F, ^ei} can be determined 
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from Eqs. ©“O- 

Table Ivni displays the blob scaling variables calculated 
with this procedure for sodium poly(styrene sulfonate) 
(NaPSS) at two different molecular weights, M = 20600 
Dalton and M = 206000 Dalton. For NaPSS, the molec¬ 
ular weight of a monomer. Mm = 206 Dalton, and the 
monomer length bm = 2.5Ai^ The intrinsic persistence 
length of NaPSS is variously estimated^^— as lying in 
the range, OA < io < IdA. For the purposes of this illus¬ 
tration, we choose a value, Iq = 12 . 5A, in order to get 
an integer number of monomers per Kuhn step. Simi¬ 
larly, a fixed value of c™ = 10 “® mol/ltr is chosen as 
the monomer concentration, since at this concentration, 
according to the phase diagram in Ref. [H, the solu¬ 
tion remains in the dilute regime for all the values of 
Nm considered here. The degree of ionization (sulfona- 
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